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Abstract 

We study the winding number dependence of the stationary states of a Bose-Einstein condensate 
in a ring-shaped lattice. The system is obtained by confining atoms in a toroidal trap with equally 
spaced radial barriers. We calculate the energy and angular momentum as functions of the winding 
number and the barrier height for two quite distinct particle numbers. In both cases we observe 
two clearly differentiated regimes. For low barriers, metastable vortex states are obtained up to a 
maximum winding number which depends on the particle number and barrier height. In this regime, 
the angular momentum and energy show, respectively, almost linear and quadratic dependences on 
the winding number. For large barrier heights, on the other hand, stationary states are obtained 
up to a maximum winding number which depends only on the number of lattice sites, whereas 
energy and angular momentum are shown to be sinusoidal functions of the winding number. 

PACS numbers: 03.75.Lm, 03.75.Hh, 03.75.Kk 
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I. INTRODUCTION 



Ultra-cold atoms in optical lattices nowadays represent a perfect laboratory system for the 
quantitative study of various condensed-matter models. In addition to offering an extremely 
high degree of experimental control and reproducibility, such systems have shown to provide 
a very useful theoretical testing ground for the quantum modeling of strongly correlated 
many-body systems. Perhaps the most simple and fecund of such applications is given 
by the dyuam.cs of dooble-we.l Bose-Einstein condensates g winch has been extensive.y 
studied in recent years |2J. In particular, for high enough barriers, a simple two- mode model 
Hamiltonian describing Josephson-oscillations and quantum self-trapping phenomena, has 
been analyzed by different authors , and later improved to better introduce interparticle 
interaction effects On the other hand, the associated dynamics has been experimentally 
observed a few years ago [8]. 

In the last years, this research work has been extended to multiple well condensates 
motivated by the possibility of generating and manipulating vortex states. In a recent 
experiment, Scherer et al. [9] have reported observations of vortex formation by merging and 
interfering multiple Bose-Einstein condensates in a confining potential. In their experiment, 
a single harmonic potential well was partitioned into three sections by barriers forming three 
independent Bose-Einstein condensates, which had no coherence between the respective 
phases. By removing the barriers, the condensates merged together and it has been shown 
that, depending on the unknown initial relative phases, the final state could acquire vorticity. 
The authors have observed single vortex states even after five seconds following the barrier 
ramp-down, indicating relatively long single vortex lifetimes. This experiment stimulated 
theoretical investigation as regards the maximum winding number a condensate in a ring- 
shaped lattice could sustain, depending on the trap properties, particularly, on the barrier 
height between wells. For example, in the case of large barrier heights, the problem of vortex 
trapping in cyclically coupled Bose-Josephson junctions has been recently addressed by 

n 

Ghosh and Sols [101 ] . They coupled initially independent Bose-Einstein condensates through 
Josephson links and allowed the system to reach a stable circulation by adding a dissipative 
term in semiclassical equations of motion. In addition, the authors analyzed the probability 
of trapping a vortex with a given winding number, finding that a necessary condition to 
obtain a metastable nonzero circulation consists of generating a winding number smaller 
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than a quarter of the total number of linked condensates. A different approach was utilized 
by Perez-Garcia et al. which by means of group-theoretical methods, studied charge 
inversion and vortex erasing in condensates confined by traps with discrete symmetries 12 ]. 
They have predicted the existence of stationary states up to a winding number which depends 
on the order of the symmetry. On the other hand, the case of a continuous symmetry, 
represented by nonrotating ring shaped traps, has attracted renewed theoretical attention 
in connection with the existence of vortex states, due to recent experimental observations of 
persistent currents in multiple connected condensates Q, Q. In particular, by computing 
the energy landscape, Capuzzi et al. [15| have shown that multiply quantized vortices and 
multivortex configurations can be locally energetically stable in a toroidal trap. 

The aim of this work is to study vortex states in a ring-shaped optical lattice with 
variable barrier heights. This will allow us to explore the transition between the limit of 
low barrier heights, where large persistent currents are expected, and the high-barrier limit, 
where Josephson-type links should rule the vortex configurations. For such purpose we 
will consider a condensate confined by superposition of a toroidal trap and equally spaced 
radial barriers forming a ring of N c linked condensates. In particular, we will investigate 
the energy and angular momentum dependence on the winding number for a wide range 
of barrier heights. Such a system could be experimentally investigated, as it corresponds 
to a relatively simple combination of previously implemented experimental settings jg. [l4]. 
Actually, new techniques have been recently developed by Henderson et al. to create 
and manipulate condensates in a variety of geometries, including the present type of ring 
lattices. 

Our work is organized as follows, in Sec. [Til we describe the theoretical framework utilized 
to generate stationary states with different winding numbers. Sec. II III is devoted to an 
analytical study of the energy and the angular momentum for limiting values of the barrier 
height. By numerically solving the Gross-Pitaevskii equation, we calculate in Sec. IIVI the 
energy and angular momentum and analyze their dependence on the winding number and 
the barrier height. Finally, a summary and concluding remarks are offered in Sec. |V] 
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II. THEORETICAL FRAMEWORK 



A. The Trap 

We consider a Bose-Einstein condensate of Rubidium atoms confined by an external trap 



V trap , consisting of a superposition of a toroidal term V tOTO and a lattice potentia 



formed by radial barriers. Similarly to the trap utilized in recent experiments 
toroidal trapping potential in cylindrical coordinates reads, 



Vl(x, y) 



13|,ll4|, the 



V toro (r,z) = - [ojy + Lo 2 z z 2 \ + V exp(-2r 2 / X 2 ) (1) 

where ui r and u z denote the radial and axial frequencies, respectively. We have set u z » u r 
to suppress excitation in the z direction. In particular, we have chosen u r /(2ir) = 7.8 Hz 
and bj z j (2tt) = 173 Hz, while for the laser beam we have set Vq = 100 %u r and Aq = 6 l r , with 



l r = yh/(muj r ). On the other hand, the lattice potential formed by N c Gaussian barriers 
located at equally spaced angular positions 9i = 2ni/N c (0 < % < N c — 1) may be written, 

w v T/ / [cos(^)i/-sin(^)x] 2 | 

V L (x } y) = V b 2^ expj 1 (2) 

for N c even. Although, for simplicity the numerical simulations of this work have been 
restricted to the case of an even number of lattice sites, most of our findings can be easily 
shown to be valid irrespective of the parity of N c . We have fixed the trap parameters at 
Xb = 0.5l r and iV c = 16. The combination of these potentials gives rise to effective barriers 
between neighboring sites, as represented in Fig. [TJ 

The Thomas- Fermi (TF) approximation, will allow us to deduce some general features of 
the condensate from the characteristics of the trapping potential. To this aim, we will take 
into account that the TF-density vanishes at points fulfilling /i < V tTa , p (x,y), \i being the 
chemical potential. On the one hand, around the trap center we can approximate Vt rap ~ 
Vo + (N c /2)Vb and, as we restrict our calculations to chemical potentials smaller than Vo, the 
condensate center will always exhibit a hole. On the other hand, between neighboring sites 
we have that the potential minimum shown in the right panel of Fig. [T], which corresponds 
to the saddle point in two dimensions, turns out to be located around x m i n = 6.6 l r , and the 
barrier height can be approximated as V min ~ ^^(6.6, 0) = 30.7 hu r + Notice that we 
have retained in this estimate only the i = term of the sum in Vl(x, y), as the exponent of 
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FIG. 1: (Color online) Left panel: equipotential curves for the external trap Vt ra , p (x, y)/fioo r with 
Vb/Tioor = 80. Right panel: effective potential barrier Vtrap(^> 0) along the positive x axis for several 
values of the parameter Vb/hu r . 

such a term vanishes, while the contributions of the remaining Gaussian terms are negligible. 
Thus, when the chemical potential becomes smaller than this potential minimum, the TF- 
density between neighboring sites vanishes. 



B. Stationary States 

Since the trapping potential is much stronger along the z axis, the motion of the particles 
in the z direction remains frozen in the ground state, and the energy per particle in the 
condensate is then given by the two-dimensional (2D) energy functional [l7] 

m = I (l^' 2 + Krap ^ + \ N9 dxdm (3) 

where ip = e l< ^|^| denotes the 2D order parameter, N the number of particles, and m the 
atom mass. The effective 2D coupling constant g = g 3 D^Jmijj z /2nh is written in terms of 
the 3D coupling constant between the atoms g^o = 47raa 2 /m, where a = 98.98 ao denotes 
the s-wave scattering length of 87 Rb, ao being the Bohr radius. In order to obtain stationary 
states, one has to perform variations of E with respect to ip keeping the number of particles 
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fixed, which yields the Gross-Pitaevskii (GP) equation 



h 2 v 2 



+ Krap + Ng \ij\ 2 V = ^ ■ 



(4) 



2m 



The ground state wavefunction is numerically obtained without imposing any constraint. 
Typical configurations are shown in Fig. [21 where we depict the equidensity curves for 
different barrier heights. The left panel corresponds to the barrier-free system, while the 
middle and right panels correspond to barrier parameters fixed at Vb/hu) r = 10 and Vb/hu r = 
80, respectively. We note that the chemical potential verifies /x ~ 77hco r > V m i n for the lower 
barrier, while the higher barrier verifies the opposite inequality fi ~ 90 hu r < V m i n . Thus, as 
we have explained in the previous section, the TF-density between neighboring sites should 
vanish in the latter configuration, while remaining finite in the former one. As regards the 
exact density, it should exhibit a similar behavior, except for a smoother decay. 

On the other hand, further constraints must be imposed in order to obtain stationary 
states with nonvanishing winding numbers, and to this aim we have employed a method 
which was described in detail in a previous work [19| . Summarizing, we begin by using the 
ground state wavefunction ip to construct an initial order parameter with vorticity [151 ]: 



where n denotes the initial winding number or topological charge. From this ansatze, it is 
straightforward to realize that ip and ipo possess the same density profile, but ip presents an 
imprinted velocity field. Such a field is irrotational everywhere, except at the origin, with 
a circulation along any closed loop around this point, which must be proportional to the 
winding number n. Next, we apply a minimizing procedure in order to obtain a stationary 
state. During the minimization such a velocity field undergoes drastic changes when large 
barrier heights are considered. Notice that the presence of radial barriers destroys the axial 
symmetry, and thus the angular momentum does no longer commute with the Hamiltonian. 
As a consequence, there is also no longer a linear relationship between the value of the 
angular momentum and the imprinted winding number. In addition, the final winding 
number n' after the minimization process may differ from the initial one n, as will be shown 
in Sec. HVl where we will also see that |n'| is bounded by a maximum value, which we shall 
call v. 




(5) 



6 




-15-) , , , , , , , , , , , | " 1 5 -| , , , , , , , , , , , 1-15-| , , , , , , , , , , , 

-15 -10 -5 5 10 15 -15 -10 -5 5 10 15 -15 -10 -5 5 10 15 

x/l r 



FIG. 2: Density isocontours for several values of the barrier parameter. From left to right: 
Vb/huJr = 0, 10, and 80. The number of particles corresponds to N = 10 5 . 

III. ANALYTICAL EXPRESSIONS 

In this Section, we shall analyze the behavior of the angular momentum and the energy 
as functions of the winding number in two limiting cases, namely in the absence of barriers 
and for large barrier heights. In both cases we will derive analytical expressions for the mean 
values of such observables in terms of ground state magnitudes. Hereafter, when referring 
to energy and angular momentum, it must be understood that they denote the mean value 
of the corresponding observables per particle. 



A. Barrier free system 

Taking into account that the ground state density exhibits a hole at the condensate 
center, for vanishing barriers and an arbitrary winding number n, the order parameter may 
be approximated as: 

Mr,0)=Mr)e M , (6) 

where we have assumed that the presence of a non-homogeneous phase in the wavefunction 
with its concomitant velocity field, does not affect appreciably the shape of the condensate 
density, as compared to that of the ground state. This occurs because the centrifugal force 
should have little effect outside the central hole, where the density is non-negligible. 
Replacing the wavefunction (Q in the first term of the integrand of ([3]) we have: 

» 2 , 2 h 2 h 2 , , 2 n 2 

_| V . W ^_ (vw2 + _ Wo)2 _ (7) 
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and one can easily obtain the mean value of the energy as a quadratic function of the winding 
number n, 



E(n) = Kn 2 + E (8) 
where E denotes the ground state energy and K is given by: 

K = j '±[Mr)} 2 rdr. (9) 

On the other hand, since the angular momentum commutes with the Hamiltonian and 
the wavefunctions ([6]) are eigenfunctions of L z , we have the linear relationship 

L z = hn. (10) 

B. High barriers 



According to Refs. [l2j and Uj the most general solution of the GP equation is a Bloch 
state of the form: 

Mr,0) = e m9 f n (r,6), (11) 

where the function f n {f,0) is invariant under a 2tt/N c rotation in 9. This type of solution 
is valid irrespective of the barrier height. The above periodicity is obvious in the modulus, 
since |/ n (r, 9)\ 2 = \ip n {T, 0)\ 2 corresponds to the periodic particle density as shown in Fig. [21 
However, the phase of ijj n (r,6) does not necessarily present this symmetry. As an example, 
in Fig. [3] we depict the phase of ^(ro,^) as a function of the angular coordinate, which 
clearly does not exhibit such a symmetry. Whereas, it is easy to verify that the phase of 
f4( r o,&) = e~* 4e ip4(r ,9) is indeed a 2tt/N c periodic function. 

In this section we will concentrate ourselves on large barrier heights. In such a limit, the 
phase turns out to be almost uniform in each site, as it may be seen from the right panel 
of Fig. [3j This is a consequence of the fact that for high barriers, the condensate portions 
become more isolated and thus the particle current diminishes. This, in turn, implies that the 
velocity field v and hence the gradient of the phase (m/h)v practically vanishes within each 
portion. Then, for large enough barriers, such that the chemical potential verifies /i < V m i n , 
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FIG. 3: Phase of the order parameter as a function of the angular coordinate for tq = 6l r , 
a winding number n = 4, and two values of the barrier height Vb/fiuj r = 60 (left panel) and 
Vb/fiUJr = 120 (right panel). The number of particles is N = 10 5 . 

the system splits into a ring of weakly linked condensates, and the order parameter can be 
approximated as follows, 

iVc-l 

Mr,0)= E y k (r,6)e i2 ™ k / N % (12) 

where \l/fc(r, 9) is a real function describing a localized state in the fc-site, i.e., having only 
a very little overlap with neighboring sites. This approximation is analogous to the tight- 
binding approximation [2l| . and the matrix element ip^Oipn of a given observable O can be 
approximated as 



N c -l „ N c -1 N c -1 



E ^kOm k +e l ~ E +e" J — E **0*fc-i, (13) 

fe=0 k=0 k=0 

where we have assumed that only the matrix elements between neighboring sites are non- 
negligible. In our notation, the site denoted by k — N c (k — —1) is identified with that of 
k = (k = N c — 1) . Then, as the summation is performed over all k values, we can rewrite 
the expression ( Fl3l) as follows: 

iVc-l 2 Nc-1 2 Nc-1 

E *fcO** +e 1 ^ E **0* fc+ i +e~ 4 ^ E (14) 

fc=0 fc=0 k=0 

Let us now assume that the observable O is the angular momentum operator: 
Using that L z is hermitian and ^/^ are real functions, the following equalities hold: 
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* fc L z ^ k = 



(16) 



V k L z ^ k+1 = -* fc+1 L 2 (17) 
Then replacing ( JIBjl and ( 1T7|) in ffT4]) we obtain 



<L^n = ~(e- - e*K) £ (18) 
and thus the mean value of the angular momentum reads: 

/N c — 1 , (j^Tj 
r n L z ip n rdrd6 = 2h sm(2nn/N c ) £ / * k —£± r rfr d6 . (19) 
fc=o 7 dB 

Note in the above expression that the winding number appears only in the argument of the 
sine, which determines that the angular momentum must exhibit a sinusoidal behavior as a 
function of n. 

In computing the mean value of the energy, it is convenient to separately treat the 
interacting and noninteracting terms. In fact, the noninteracting part, which corresponds to 
the observable O s = — ^r~ + Vt rap , may be treated analogously to the angular momentum, 
yielding 

^ k O s ^ k+1 = V k+1 Os V k , (20) 

which replaced in fTH|) leads to 



N c -1 



f I 'il>*d,il> n rdrd6= £ / fy k d s V k rdrdd 

N c -1 

+ 2 cos(27m/iV c ) J2 / ^kOs ^k+i rdrdO. (21) 

fc=0 

On the other hand, the interacting term corresponds to the nonlinear operator O n i = 
7;N\tp n \ 2 , proportional to the particle density, for which the approximation (|14p is no longer 
valid. In this case, the mean value must be computed to first order in the product of adjacent 
wavefunctions ^f^k+ii yielding 

-N f f \ip n \ 4 rdrd6= Y, I I -N^irdrdO 
2 J J k=Q J J 2 
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N c -1 



+ 2gNcos(27rn/N c ) £ I /(*k*k+i + *fc^+i) r drdQ. (22) 

fc=0 

Finally, from (I2ip and ( 122|) we may conclude that the mean value of the energy will present 
a dependence of the type cos(27rn/iV c ) on the winding number. 

IV. ENERGY AND ANGULAR MOMENTUM: NUMERICAL RESULTS 
A. Low barriers 

A vanishing barrier (V& = 0) yields a ground state chemical potential fi = 33.19 hu r 
(fi = 73.31 hjj r ) for iV = 10 3 (N = 10 5 ) particles. Such values turn out to be smaller than 
Vb, and thus, as discussed in Section II, the correspondin g co ndensates exhibit a hole at 



their centers 20(. It has been shown, both experimentally 13( and theoretically 15J, that 
this type of configuration forming a torus can host locally stable multiply quantized vortices. 
We have found such stable vortices up to a maximum winding number of v = 7 ( v = 9 ) for 
iV = 10 3 (N = 10 5 ), reflecting the fact that v depends on the number of particles. Actually, 
we will next see that v also depends on the barrier height. 

In the top panel of Figs. HJand |5]we depict the excitation energy E' n = E{n) — E as a 
function of the winding number. In the case of a vanishing barrier, we have also numerically 
evaluated the parameter K of Eq. (Q, obtaining the values 0.01156 and 0.01114 for iV = 10 3 
and iV = 10 5 , respectively. Thus, we have drawn the curve arising from formulae OH]), which 
shows an excellent agreement with the numerically obtained values (circles), as seen from 
Figs. S] and El It may also be observed that the energy decreases for increasing values of 
the barrier height, retaining, however, the approximately quadratic behavior predicted by 
Eq. OH]). 

The bottom panel of Figs. H] and [5] shows the angular momentum as a function of the 
winding number. Note that for a vanishing barrier, the numerical results correspond to 
the straight line predicted by Eq. f llOp . while we may observe that the angular momentum 
decreases and curves down for increasing values of the barrier height. In the same sequence, 
we may also observe a decreasing value of the maximum winding number v the system 
can host. In fact, we have found that for a given barrier height, the time evolution with 
a moderate dissipation of an initial state with winding number n > u, may lead to 
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FIG. 4: Excitation energy per particle E' n (top panel), and angular momentum per particle 
L z = (L z ) (bottom panel), as functions of the winding number for N = 10 3 particles and Vb/huj r = 
(circles), 3 (stars), and 5 (triangles). The solid lines correspond to Eq. fl5J) (top) and Eq. (TH)|) 
(bottom). 

different processes of vorticity reduction, which involve vortex motion along the barriers. 
Such a vortex dynamics may either simply consist of a vortex escaping phenomenon or, by 
contrast, it may contain more complex processes such as, vortex escaping with the aid of 
antivortices. We describe in the following this last particular process. Fig. E] shows a series 
of snapshots in time of the phase distribution when a winding number n = 9 is initially 
imprinted in the condensate of 10 5 particles and barrier height Vb/hu r = 40 (panel(a)). We 
recall that the maximum winding number for this barrier height is v — 7, as shown in Fig. |5j 
The time evolution first yields a ring of 16 vortices of topological charge m — +1, denoted as 
white dots in panel (b), which move along the barriers from the condensate center outwards. 
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FIG. 5: Same as Fig. U]for the condensate of N = 10 5 particles with barrier heights Vb/hco r = 
(circles), 10 (squares), 40 (stars), and 60 (triangles). 

Thus, the circulation along the dashed line circle in panel (b) still bears n = 9, while the 
same calculation for a circle within the vortex ring would yield -7. However, to reach the 
condensate border, each vortex should overcome the high density region corresponding to 
the potential minimum shown at the left panel of Fig. [TJ which approximately lies around the 
dashed line circles shown in Fig. |6j Such a displacement would produce an overall increase in 
energy that is actually avoided by the generation of an outer ring of antivortices, denoted as 
black dots in panel (c). This leads to the formation of vortex-antivortex pairs, i.e., dipoles, 
of lower kinetic energy that are then able to keep on moving towards the condensate edge, 
leaving behind them a pattern of centered vorticity of winding number -7. This process is 
accompanied by a progressive reduction of the separation of both charges at each dipole. In 
the last panel (d), we may observe that each dipole is about to annihilate itself, while the 
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FIG. 6: (Color online) Snapshots of the phase distribution in the condensate of 10 5 particles with 
a barrier height Vb/huj r = 40 and an initially imprinted winding number n = 9. Panels (a) to (d) 
correspond to uj r t = 0, 0.374, 0.796, and 3.183, respectively. White dots in panels (b) to (d) denote 
vortices with topological charge m = +1, while black dots in (c) and (d) denote antivortices with 
m = — 1. The dashed line circles in (b) and (d) correspond to circulations yielding initial n = 9 
and final n' = —7 winding numbers, respectively. 

final n! = —7 central vorticity pattern becomes evident by following the circulation along 
the dashed line circle. 

Finally, we want to remark that time evolutions with n < v always lead to an identical 
final winding number n' = n, with the vorticity concentrated at the central hole. 
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FIG. 7: Excitation energy per particle E' n (top panel) and angular momentum per particle L z = 
(L z ) (bottom panel), in units of their maximum values E max and L max , respectively, as functions 
of the winding number for the condensate of N = 10 3 particles, with barrier heights Vb/huj r = 10 
(crosses) and 40 (squares). The solid lines represent the functions: [1 — cos(27rn/iV c )]/2 (top panel) 
and sm(2irn/N c ) (bottom panel). 

B. High barriers 

Now we will study the numerical results obtained for the high-barrier regime which was 
analytically investigated in Sec. IIII Bl We depict in Figs. [7] and El the excitation energy 
(top panel) and the angular momentum (bottom panel), as functions of the winding number 
for the condensates of 10 3 and 10 5 particles, respectively. We may observe that while both 
magnitudes are well described by the periodic functions obtained in IIII B\ there still exist 
for the lower barrier heights, Vb/hu r = 10 and 80, some noticeable differences between the 
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FIG. 8: Same as Fig. [7] for the condensate of N = 10 5 particles with barrier heights Vb/hoj r = 80 
(crosses) and 120 (squares). 

numerical values and the predicted curves. 

Irrespective of the number of particles, we have observed in this regime the processes 



predicted by Perez-Garcia et al. 11], namely the so-called charge erasure, which occurs 
only for an even number of lattice sites, and charge inversion. In fact, for an imprinted 
winding number n = N c /2 = 8, the phases between neighboring sites differ in tc (top panels 
of Fig. |9]), leading to a vanishing net flux of particles, or in other words, to an erasure of 
vorticity. Notice that such a state presents a vanishing angular momentum and a maximum 
of energy, as shown in Figs. [7] and [HJ On the other hand, charge inversion takes place if 
n > N c /2, as shown in the bottom panels of Fig. [9j In this case, the phase difference between 
neighboring sites A0 becomes larger than tt, which amounts to an effective phase difference 
of 27T — A(p. Thus, the phase derivative dip/ 89 changes sign, giving rise to an inverted 
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current density with a corresponding negative angular momentum (see Figs. [7] and [8]), and 
a final winding number nf = n — N c . We have found that repeating the dissipative real-time 
simulation of Sec. IIV A\ but for a higher barrier of Vb/hu r = 80, leads to vortex escaping 
along each of the 16 barriers and a final state of vorticity n' = —7, localized within the 
condensate central hole. However, it is interesting to notice that in contrast to the low 
barrier regime (see Fig. [6]), where the vortices may be viewed as point charges, we observe 
in Fig. [10] a remarkable phenomenon of vorticity spreading of the escaping vortices. 

Summarizing, we may realize that there exists in this regime, irrespective of the number 
of particles, N c = 16 different stationary states, namely N c — 2 = 14 'multivortex' states 
corresponding to final winding numbers 1 < \n'\ < N c /2 — 1 = 7, plus a single excited 
state of maximum energy and vanishing particle current corresponding to n — N c /2 = 8, 
and the ground state ( n — 0). These results are valid for any even number of lattice sites. 
For an odd number of sites, on the other hand, the only difference would be given by the 
absence of the single excited state of maximum energy and vanishing angular momentum, 
i.e., we would have N c — 1 'multivortex' states correspondin g to final winding numbers 
1 < \n'\ < {N c — l)/2, apart from the ground state with n = Hi, \12 \. 

Finally, in Fig. [TTJwe depict the chemical potential of the ground state as a function of the 
barrier height Vj,, jointly with the effective potential barrier minima V m i n shown in Fig. [TJ 
We notice that V min approximately intersects the chemical potentials corresponding to 10 3 
and 10 5 particles at Vt, = 5 Tn>j r and = 60 hu r , respectively. Such intersections correspond 
to transitions between the low- and high-barrier regimes above described, i.e., they indicate 
the barrier height above which the quantum tunneling regime becomes dominant. 



V. SUMMARY AND CONCLUDING REMARKS 

We have analyzed the angular momentum and energy as functions of the winding number 
in a ring-shaped lattice, finding two clearly differentiated regimes, according to the relative 
value of the barrier height as compared to the chemical potential. On the one hand, we have 
found that for low barriers, there exists an almost quadratic dependence on the winding 
number for the energy and a linear dependence for the angular momentum, as it is ana- 
lytically derived in the absence of barriers. For high barriers, on the other hand, we have 
found that the energy and the angular momentum turn out to be accurately described by 
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FIG. 9: (Color online) Phase of the order parameter as a function of the angular coordinate (both 
in units of it) for tq = 6l r (left panels), and phase distribution in the x-y plane (right panels), in 
the condensate of 10 5 particles with a barrier height Vb/huj r = 80. The imprinted winding number 
is n = 8 (top panels) and n = 9 (bottom panels). The long (blue) vertical arrow in the left-bottom 
panel shows a phase difference between neighboring sites larger than n, which amounts to the 
smaller phase difference denoted by the short (red) arrow. Note also the slightly negative phase 
slope within each site, d(p/d9 < 0, which yields the final winding number n' = 9 — N c = —7 and 
the corresponding negative angular momentum, as seen in the bottom panel of Fig. 

sinusoidal functions, as predicted from a simple analytical model. We have observed these 
regimes irrespective of the number of particles, and a remarkable difference between them 
consists in that for the low-barrier regime, metastable vortex states are obtained up to a 
maximum winding number u, which depends on the particle number and the barrier height, 
whereas, for high barriers, such states are obtained up to a maximum winding number which 
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FIG. 10: (Color online) Snapshot of the phase distribution in the condensate of 10 5 particles with 
a barrier height of Vb/huir = 80 and an initially imprinted winding number n = 9. Real-time 
dissipative evolution at t = 0.2387 u~ showing vorticities with a topological charge m = +1 
spread along each barrier. Such a spreading may be appreciated by following the circulation along 
the rectangular contour enclosing one of the barriers. 

only depends on the number of lattice sites. Another difference between both regimes stems 
from the way in which vorticity becomes reduced from an initially imprinted winding num- 
ber exceeding v. For a low barrier, this process involves the displacement outwards along 
the barriers of positive point charges, jointly with an eventual generation of negative charges 
at the condensate edge. Such antivortices penetrate along the barriers to form dipoles that, 
by diminishing the energy from velocity field, are then able to cross the maximum density 
region outwards. On the other hand, the vorticity reduction in the large-barrier regime 
proceeds through the appearance of stretched positive vorticities along each barrier, which 
eventually escape from the condensate. 

To conclude, we believe that the present study could also be regarded as a first stage 
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FIG. 11: Ground state chemical potential /i of the condensate containing 10 5 particles (squares) 
and 10 3 particles (triangles) as a function of the barrier height Vj,. The solid line corresponds to 
the effective potential barrier minima of Fig. [TJ All quantities are given in units of %uj t . 

to establish a Bose-Hubbard model for ring-shaped optical lattices with a large occupation 
number per site. 

Acknowledgments 

DMJ acknowledges CONICET for financial support under Grant No. PIP 
11420090100243. 



[1] R. Gati and M. K. Oberthaler, J. Phys. B 40 R61 (2007). 

[2] M. Asad-uz-Zaman and D. Blume, Phys. Rev. A 80, 053622 (2009); N. Teichmann, M. Es- 

mann, and C. Weiss, Phys. Rev. A 79, 063620 (2009); D. Witthaut, F. Trimborn, and S. 

Wimberger, Phys. Rev. A 79, 033621 (2009); M. T. Martinez, A. Posazhennikova, and J. 

Kroha, Phys. Rev. Lett. 103, 105302 (2009); B. Julia-Diaz, D. Dagnino, M. Lewenstein, J. 

Martorell, and A. Polls, Rev. A 81,023615 (2010). 
[3] I. Zapata, F. Sols, and A. J. Leggett, Phys. Rev. A 57, 28 (1998). 
[4] I. Zapata, F. Sols, and A. J. Leggett, Phys. Rev. A 67, 021603 (2003). 
[5] A. Smerzi, S. Fantoni, S. Giovanazzi, and S. R. Shenoy, Phys. Rev. Lett. 79, 4950 (1997). 



20 



[6] S. Raghavan, A. Smerzi, S. Fantoni, and S. R. Shenoy, Phys. Rev. A 59, 620 (1999). 
[7] D. Ananikian and T. Bergeman, Phys. Rev. A 73, 013604 (2006). 

[8] M. Albiez, R. Gati, J. Foiling, S. Hunsmann, M. Cristiani, and M. K. Oberthaler, Phys. Rev. 
Lett. 95, 010402 (2005). 

[9] D. R. Scherer, C. N. Weiler, T. W. Neely, and B. P. Anderson, Phys. Rev. Lett. 98, 110402 
(2007). 

[10] P. Ghosh and F. Sols, Phys. Rev. A 77, 033609 (2008). 

[11] V. M. Perez-Garci'a, M. A. Garcfa-March, and A. Ferrando, Phys. Rev. A 75, 033618 (2007). 
[12] A. Ferrando, Phys. Rev. E 72, 036612 (2005). 

[13] C. Ryu, M. F. Andersen, P. Clade, Vasant Natarajan, K. Helmerson, and W. D. Phillips, 

Phys. Rev. Lett. 99, 260401 (2007). 
[14] C. N. Weiler, T. W. Neely, D. R. Scherer, A. S. Bradley, M. J. Davis, and B. P. Anderson, 

Nature 455, 948 (2008). 
[15] P. Capuzzi and D. M. Jezek, J. Phys. B 42, 145301 (2009). 

[16] K. Henderson, C. Ryu, C. MacCormick, and M. G. Boshier, New J. Phys. 11, 043030 (2009). 
[17] Y. Castin and R. Dum, Eur. Phys. J. D 7, 399 (1999). 

[18] E. P. Gross, Nuovo Cimento 20, 454 (1961); L. P. Pitaevskii, Zh. Eksp. Teor. Fiz. 40, 646 

(1961) [Sov. Phys. JETP 13, 451 (1961)]. 
[19] D. M. Jezek, P. Capuzzi, and H. M. Cataldo, J. Phys. B 41, 045304 (2008). 
[20] H. M. Cataldo and D. M. Jezek, Eur. Phys. J. D 54, 585 (2009). 

[21] N. W. Ashcroft and N. D. Mermin, Solid State Physics, (Saunders College Publishing, Fort 
Worth, 1976), chap. 10. 

[22] We have employed a phenomenological damped GP equation (see, e.g., K. Kasamatsu, M. 
Tsubota, and M. Ueda, Phys. Rev. A 67, 033610 (2003)) with a damping coefficient 7 = 0.03. 



21 



